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Abstract 

We consider a discrete model of euclidean quantum gravity in four dimensions based 
on a summation over random simplicial manifolds. The action used is the Einstein- 
Hilbert action plus an i? 2 -term. The phase diagram as a function of the bare coupling 
constants is studied in the search for a sensible continuum limit. For small values 
of the coupling constant of the R 2 term the model seems to belong to the same 
universality class as the model with pure Einstein-Hilbert action and exhibits the 
same phase transition. The order of the transition may be second or higher. The 
average curvature is positive at the phase transition, which makes it difficult to 
understand the possible scaling relations of the model. 



Partly supported by the KBN grant no. 2 0053 91 01 



1 



1 Introduction 



Understanding the theory of quantum gravity remains one of the greatest challenges 
in theoretical physics. One can try to circumvent the problem by embedding gravity 
in a larger theory like string theory. This is in many respects an appealing approach, 
but it seems to have lost some of its momentum and there is not much hope that it 
will be possible in the near future to arrive in a natural way to an effective theory 
of gravity (and matter) in four dimensions. Discussing basic principles one should 
maybe not be so worried about our technical inability to deduce the consequences 
of string theory since this is not the first time in theoretical physics we are unable to 
extract anything but the simplest perturbative consequences of an otherwise healthy 
and probably correct theory. It is more worrisome that string theory is not (yet) well 
defined beyond the loop expansion. Again this might not seem so disastrous since 
the same was (and to some extent still is) true for ordinary field theory. However, 
in the last three years there has been a significant progress in our understanding 
of the problems connected with the summation of all loops in string theory. The 
message has not been encouraging. At the moment we have no general principles 
which allow us to define in an unambiguous way the summation over all genera in 
string theory. From this point of view it might be somewhat premature to announce 
string theory as the fundamental theory. Strictly speaking it is not yet a theory but 
a set of rules which allow us to calculate certain perturbative quantities. 

If we decide to drop string theory as the theory which will teach us the nature of 
quantum gravity it might be (good ?) conservative policy to stay entirely within the 
ordinary field theoretical framework. At first glance it does not look too promising 
either. The four dimensional theory is hampered by being non-renormalizable and 
we do not at present know any example where such a theory can be defined non- 
perturbatively, is non-trivial and at the same time satisfies what we usually view 
as the axioms of quantum field theory. At the same time one of the lessons of the 
last thirty years is that field theory is deeply connected to the theory of critical 
phenomena via the path integral and has a natural formulation in Euclidean space- 
time. But precisely for gravity the Einstein-Hilbert action is unbounded from below 
due to the conformal mode and the Euclidean path integral is ill-defined. One could 
try to make sense of the unboundedness of the action either by a contour rotation 
associated with the conformal mode as suggested by Hawking and others or by 
stochastic regularization (the so-called fifth time action) as advocated by Greensite 
and Halpern 0, [§. It is not the purpose here to enter into a discussion of the virtues 
and drawbacks of these interesting suggestions. Let us only mention that they do 
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not really have the flavour of general descriptions based on ordinary field theory. Of 
course it is wise to bear in mind that if any field theory should depart from basic 
axiomatic principles it is quantum gravity, but lacking a general alternative we have 
decided to return to analysis of quantum gravity in the context of ordinary field 
theory. 

Field theory suggests one rather simple minded way out of the above mentioned 
problems. This was already discussed long ago by Weinberg who called it asymptotic 
safety. The idea is simply that when we, by means of the renormalization group 
equations, work our way back from the infrared fixed point where the Einstein- 
Hilbert action seems a good effective description we will at some point reach a non- 
trivial ultraviolet fixed point. In addition the associated critical surface is assumed 
to be finite dimensional, which means that only a finite number of parameters are left 
arbitrary in the theory, which from this point of view can be said not to differ much 
from ordinary renormalizable field theories. The effective Lagrangian description of 
the theory by means of fields suitable for the infrared fixed point might then be an 
infinite series 



A- j^ R + h R2 + f^RpuR^ 



which might even be non-polynomial, but which might now (and we will assume this 
is the case) make sense if we make a formal rotation to Riemannian space where 
the metric has the signature {+1,+1,-|-1,-|-1} (the generalization of the rotation to 
Euclidean space in ordinary field theory). 

Of course one weakness in this scenario is that the existence of the ultraviolet 
fixed point has been entirely hypothetical. Further we have not exactly been flooded 
with examples of this kind in field theory, as already mentioned. Finally, and we 
agree that this point is an annoyance, such a solution does not have the appeal and 
aesthetical beauty of the original theory which we want to quantize. If the quantum 
theory of gravity offers to us a solution like ( |1.1| ) the next task must be to find a 
simpler description in terms of other variables, maybe somewhat like the switch in 
hadronic physics from hadrons to quarks and gluons. 

A few things have happened since Weinberg outlined the above strategy. There 
exist now regularizations of the path integral which allow us to define theories like 
( |1.1|) non-perturbatively and further the progress in computer science has made it 
possible to calculate approximately these path integrals. It is therefore possible to 
explore, by numerical means, the phase diagram of the regularized theory and try 
to locate phase transition points in the coupling constant space. If the transitions 
are of second order one can attempt to define a continuum limit. The approach has 
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one virtue: it requires only a finite amount of work to verify whether the idea of a 
non-trivial ultraviolet fixed point is viable or not. 

The rest of this article is organized as follows: In section 2 we define the dis- 
cretized model to be used to regularize the path integral. Section 3 provides some 
details about the numerical method used, while section 4 describes the results of ex- 
tensive numerical simulations. Section 5 contains a discussion of the results obtained 
so far. 

2 The discretized model 

The continuum theory of gravity is reparametrization invariant. If we discretize the 
theory in order to regularize it we will have to break this invariance provided the 
action depends on the metric. An alternative is to consider theories which depend 
only on topology. A very interesting approach in this direction in three dimensions 
has recently been followed by many people following the work of Turaev and Viro || 
and it has now been generalized to four dimensions ||. Unfortunately the precise 
connection with the usual continuum version of Einstein gravity is not yet clear, 
especially in four dimensions. 

If we restrict ourselves to the conservative approach of discretizing Einstein's 
theory of gravity we will break reparametrization invariance since an action like 
fll.ip depends explicitly on the metric and without having done anything yet we can 
already now say that the most important question to answer in case one manages 
to carry out successfully the program outlined in the introduction will be to check 
that the theory defined by approaching in a well defined way the ultraviolet fixed 
point is really reparametrization invariant. This is by no means obvious since the 
regularization has broken this invariance explicitly. 

2.1 Quantization of Regge calculus 

Two rather different regularization schemes have been suggested. The oldest one 
goes back to Regge and we will call it Regge calculus. It was originally invented 
as a means to approximate a given smooth manifold by a piecewise flat Rieman- 
nian manifold, obtained by a suitable triangulation of the smooth manifold. Regge 
showed that it was still possible to assign in a sensible way a concept of curvature 
to such a piecewise flat manifold. For d- dimensional manifolds the building blocks 
would be d- dimensional simplices and the curvature assigned to the d— 2- dimensional 
sub-simplices. In this way one has both volume and curvature assigned to the piece- 
wise flat manifold and it is possible to approximate the continuum Einstein-Hilbert 
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action which reads 

S[g] = \J d d i^g - ^ J d d ^R (2.1) 
by the following discretized expressions: 

Jfty/g - H V M) (2-2) 

3 

[d d ^R - J2 V M 

J j 

In ( |2.2j )-( |2~3| ) the summation is over all d — 2-dimensional sub-simplices j with vol- 
ume Vj(d — 2). <5j is the so-called deficit angle associated with the d — 2 dimensional 
sub-simplex j, while the ^-dimensional volume Vj(d) associated with the d — 2 di- 
mensional sub-simplex j is defined as 

W<0 = 5(3TTjgW> P-4> 

where the summation is over all <i-dimensional simplices i which contain the subsim- 
plex j. It is possible to show that by a suitable refinement of the triangulation of the 
smooth manifold the discretized expressions (|2.2j ) and ( [2. 3D will actually converge to 
the continuum value. In this way Regge calculus provides a geometrical coordinate- 
independent description of gravity where it is natural to use the length of the links 
(the geodesic length of the edges in the triangulation of the given manifold) as dy- 
namical variables since they completely specify the flat (i-simplices used as building 
blocks. Originally the method was used mainly in a classical context where there 
are no conceptional problems connected with the approach. However, already as 
early as in 1968 Regge and Ponzano in an impressive paper, which contains also 
the seed to recent development in topological gravity mentioned above, pointed out 
that quantum mechanical amplitudes in three-dimensional Regge calculus can be 
defined by a functional integral and maybe computed non-perturbatively. However, 
if we seriously want to apply the Regge calculus directly in the functional integral it 
looses some of its beauty. In the classical context the geometry was specified by the 
length of the links of the building blocks and the incidence matrix which specified 
how the building blocks were glued together. This incidence matrix which deter- 
mines the topology was fixed and not considered a dynamical variable. When we use 
the Regge formalism in the path integral the situation is opposite. The task is not 
to approximate a given continuum Riemannian manifold but (at least) to sum over 
equivalence classes of metrics associated with a given manifold. Unfortunately there 
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Vj(d 



(2.3) 



5 



is no simple one-to-one correspondence between link length and equivalence classes 
of metrics, as is easily seen by considering triangulations of the two-dimensional 
plane. Obviously many triangulations correspond to the same Riemannian geome- 
try. This means that a nasty jacobian is involved if we want to use link lenghts as our 
integration variables. In addition one has to choose an integration measure which 
ensures that the link lengths satisfy the triangle inequalities and their higher di- 
mensional analogues, which express that the fc-dimensional volume of /c-dimensional 
sub-simplices in the given triangulation must be positive. A great deal of work has 
gone into understanding and repairing these shortcomings of conventional Regge 
calculus. For a recent excellent review and references we refer to ||. While the 
classical Regge calculus gives a coordinate independent geometrical description of 
gravity it of course has nothing like reparametrization invarianceQ. It is therefore 
necessary to prove that this invariance is recovered at the point in coupling con- 
stant space where the continuum limit is taken. Unfortunately the "quantum Regge 
calculus" has not, in our opinion, been so successful in this respect. For instance 
computer simulations in [|llj seemingly give the wrong coupling to Ising spins in two 



dimensions where the coupled spin-gravity system can be solved explicitly in the 
continuum. Hopefully this is due to problems with the simulations rather than a 
basic flaw in the approach, but we do not know for sure. 

2.2 Dynamical triangulation gravity 

Due to the above mentioned problems with the translation of classical Regge calcu- 
lus to quantum theory we will here use another approach which has been extensively 
used in the last few years in the study of two-dimensional gravity and non-critical 
strings (which are nothing but two-dimensional gravity coupled to special matter 



fields) |]12| , y| [14], jlSj . We will call it dynamically triangulated gravity or (inter- 
changeably) simplicial gravity. 

In this approach the fundamental building blocks are regular simplices. In d = 2 
this means equilateral triangles, in d = 3 regular tetrahedra. One now constructs the 
manifolds by gluing together the regular <i-dimensional simplices along their d — 1 
dimensional sub-simplices, in such a way that they form a piecewise flat manifold. 
The dynamics is shifted from the length of the links to the connectivity of the 
piecewise linear manifold and as we shall see there will not be the over-counting 
present in the quantum Regge prescription. 

The assignment of volume and curvature for a given triangulation T created by 



2 One might try to define the analogue of local coordinate transformations, see for instance 
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and 



gluing together the regular simplices is in this case very simple. Let us introduce the 
following notation: An i-dimensional (sub)-simplex is denoted rii, i.e. vertices are 
denoted no, links n\ etc. The total number of such (sub)-simplices in T is denoted 
Ni(T). By the order o(rii) we understand the number of d-dimensional simplices 
which share the sub-simplex nj. We will usually consider only the class of regular 
simplicial manifolds where we have put the following restrictions on o(nj): 

o(ri d _ 1 ) = 2, o(ni)>d-i + l, (i<d-2). (2.5) 

- 12.!) reduce to 

J d d £^g oc N d (2.6) 

/ d d ^Rcx E(cd-o(n d _ 2 )). (2.7) 

The constant c d in Q2.7| ) should be adjusted in such a way that for a hypothetical 
triangulation of flat space the sum should give zero. For d = 2 one can triangulate 
flat space with regular triangles. The order of each vertex is 6 and consequently 
C2 = 6. Higher dimensional flat space does not admit a regular tessellation, but one 
can still ask for the average value of o(n d _2) required to fill up ci-dimensional flat 
space. The angle 9 d between two d—1 dimensional simplices belonging to the same 
(i-simplex is given by 

cos9 d = l/d (2.8) 
and in order to fill up <i-dimensional space we need to have 

o(n d _ 2 ) = 2n/9 d = c d (2.9) 

This is the constant which enters in (|2.7|) . We find 

c 2 = 6, c 3 = 5.104, c 4 = 4.767. (2.10) 

Let us further note that 

E oM = { ^p±N d (2.11) 

since there are ^ \ J d—2 dimensional sub-simplices in a <i-dimensional simplex. 

The discretized version of the continuum action can now, for a given triangulation 
T, be written as 

S d [T] = K d N d (T) - K d . 2 N d _ 2 (T). (2.12) 
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A more general action would be the following: 

d 

S d {T} = J2K l N l (T) (2.13) 

involving the fugacities for all different z-dimensional sub-simplices. Of course one 
can choose to consider actions which cannot be expressed entirely in terms of the 
iVj's. The higher derivative terms which can be added to the continuum action ( |2.1| ) 
will in general be of this kind, and we are going to consider them later, but let 
us for the moment discard such terms. Not all iVj's are independent. The relations 
between the iVj's can be worked out by the requirement that the triangulation should 
be locally homeomorphic to R d . This means for instance that all n^'s sharing a given 
vertex no should be homeomorphic to the unit ball in R d . Similar restrictions hold for 
the neighbours to an i-dimensional simplex n« in the triangulation, and the relations 



this imposes on the N^s are summarized in the so-called Dehn-Sommerville fl6 
relations 

^ = £(-l) fe+d ( - + W ( 2 - 14 ) 

k=i \ I 

valid for all % > 0. These relations are not independent, but allow us to eliminate 
all A^j+i's if d is even and all A^'s if d is odd. In the case of even dimensions we 
have for a given triangulation T in addition Euler's relation 

it(-iy +d N t (T)=x d (T). (2.15) 

i=0 

where Xd{T) denotes the Euler characteristic of the piecewise linear manifold which 
corresponds to the triangulation T. Of course this relation is only useful if we know 
the Euler characteristic of the triangulation T. As we shall see the restriction of 
topology is very important, and in case we fix the topology of T we can use ( |2.15|) 
to eliminate for instance N . In case the topology is not fixed we can trade N (T) 



for Xd(T). For odd dimensions (|2.15|) follows from the Dehn-Sommerville relations 
with 

Xd=2n + 1=0. (2.16) 

In odd dimensions the Euler characteristic is identically zero for any simplicial man- 
ifold, and it follows just from the requirement of local homeomorphism to R d . 

The recipe for going from the continuum functional integral to the discretized 
one is now: 



,»■, * £ (2.17) 

J T 1 . /T 
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The formal integrations on the lhs of ( |2.17|) and ( |2.18|) are over all equivalence classes 
of metrics, i.e the volume of the diffeomorphism group is divided out. T denotes 
a suitable class of triangulations. One class of constraints is given by Q2.5p, but 
it should be stressed that such short distance restrictions are not expected to be 
important in the scaling limit. 

Since different triangulations give rise to different curvature assignment one can 
view the above summation as a summation over different Riemannian manifolds. 
There is no problem with over-counting in this formulation. The idea of the contin- 
uum functional integral is precisely to perform such a sum with weight e~ s ^. Of 
course the discretized sum on the rhs of ( [2.17| ) and ( |2.18| ) can only be viewed as 



an approximation to the continuum expression which hopefully "converges" in the 
scaling limit to the correct expression. The questions which are difficult to answer 
are whether the class of piecewise flat manifolds is "close" to the class of Rieman- 
nian manifolds and whether the piecewise linar manifolds are selected sufficiently 
uniformly with respect to Riemannian manifolds that Q2.17D and ( |2.18| ) are good 
approximations. Unfortunately there is no weak coupling expansion where one can 
check this, but it is very encouraging that the formalism is known to work in the 
two-dimensional case, even if one couples conformal matter with central charge c < 1 
to the system. In this case it is possible to solve both the continuum and the dis- 
cretized system. In particular we see that reparametrization invariance is recovered 
in the scaling limit. 

Since d = 4 has our main interest in this work we can write our partition function 

as 

Z(k2,Ka) = J2 e~ KiNi+K2N2 (2.19) 
TeT 

This is the grand canonical partition function, where the volume of the universe 
can vary. It is sometimes convenient to change from the grand canonical ensemble 
described by ( |2.19| ) to the canonical ensemble where the volume iV 4 is kept fixed. 
The corresponding partition function will be 

Z(k 2 ,N a ) = ]T e K27V2(T) (2.20) 
T~r(iv 4 ) 

Z(k 2 ,k 4 ) = £ Z(K 2 ,N 4 )e- KiN \ (2.21) 

N 4 

If the entropy, i.e. the number of configurations for a given iV 4 , is exponentially 
bounded it is easy to prove that there is a critical line k 4 = k c a (k 2 ) in the (k 2 , k 4 )- 
coupling constant plane. For a given k 2 the partition function ( |2.19|) will then be well 
defined for k 4 > k\{k 2 )- Let us call this domain in the coupling constant plane T>. 
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Critical behaviour can be found only when we approach the boundary dT> which we 
denote the critical line, but in general we only expect interesting critical behaviour 
at certain critical points on &D (i.e. at certain values of k 2 ). These are the points 
we are looking for in the numerical simulations. 

Let us end this subsection by discussing a point which is worth emphasizing. 
We have been deliberately vague defining the class of triangulations T(A r 4 ) over 
which the summation is to be performed in a formula like (|2.20|) . Already in two 
dimensions where the classification of topology is so simple (it is defined by the 
Euler number \) an unrestricted summation over manifolds of different topology 
does cause problems. In fact the two-dimensional analogue of (|2.20|) and (|2.21|) 
which by means of Euler's relation can be written: 

Z{k ,N 2 ) = Yl e ' K ° X ( 2 - 22 ) 

T~T(7V 2 ) 

Z(R ,R 2 ) = Y, Z ( k o,N 2 )e- k ^ (2.23) 

does not make any sense. The well known reason is that the number of trian- 
gulations, Af(N 2 ), which one can make by gluing together a given number N 2 of 
equilateral triangles to a two-dimensional surface is too large. It grows factorially 
fast: A/"(AT 2 ) > N 2 \. This means that the two-dimensional analogue of fl2.21| ) is never 
convergent. This is not a spurious result of a perverse discretization. An analogous 
result has been proven in the continuum two-dimensional theory where the volume 
of moduli space grows at least factorially with the genus |17|]. It is the same effect 
we observe in the discretized case: As long as we fix the topology, i.e. the Euler 
characteristic \ °f the surface, we have: J\f x (N 2 ) ~ N 2 X 3 exp(k 2c N 2 ), i.e. only an 
exponential growth of the number of surfaces. In this case ( |2.23| ) is well defined 
for a certain range of k 2 's. However, an unrestricted summation over topologies 
makes the subleading pre-exponential factor N 2 X dominant when \x\ ~ N 2 since 
7x 5 x/4. 

In higher dimensions the situation is of course only worse and the best we can 
hope for is a well defined expression for a fixed (or at least restricted) topology. 
In the following we will restrict ourselves to four-dimensional manifolds with the 
topology of S 4 . 

The above outlined non-perturbative definition of gravity has nothing to add to 
our understanding (or rather lack or understanding) of the question of whether or 
not to sum over topologies in quantum gravity. Apart from the problem that the 
topologies of non-simply connected four- dimensional manifolds cannot be classified 
in a sensible way, the partitionfunction does not even make sense if we only restrict 
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ourselves to the sub-class of topologies which one can construct by simple analogy 
to two-dimensional surfaces of genus g. In the rest of this article we will only be 
interested in the search for non-trivial fixed point of the above defined theory (mod- 
ified with higher curvature terms) where the class of triangulations T corresponds 
to manifolds with the topology of S 4 . 



2.3 Higher curvature terms 

There is no straightforward generalization of Regge's work to theories of gravity 
which involve higher derivative terms like R 2 in the action. The reason is that 
Regge viewed the piecewise flat manifold, not as a discrete approximation to an 
underlying continuum surface, but as one where curvature could be defined in a 
mathematical stringent way, entirely in terms of the geometrical concepts involved in 
parallel transportation. In this way the curvature occurs in 5-functions on the lattice 
geometry, with support on the d — 2-dimensional sub-simplices. Anything other than 
the Einstein action (and the cosmological term) will then involve higher powers of 5- 
functions and this means that for piecewise flat manifolds, interpreted as by Regge, 
terms like / ^JgR 2 are infinite. In order to make sense of higher derivative terms one 
has to change the perspective on Regge calculus somewhat as advocated by Hamber 



and Williams [18| and view the lattice geometry as representing an approximation 
to some smooth geometry and the local curvature as some average curvature for 
a small volume. In fact our formulas for Regge calculus have already hinted this 
interpretation in the sense that we have assigned a volume density V nd _ 2 (d) to each 
d — 2-dimensional sub-simplex 72^-2 which can be viewed as an appropriate share of 
the volumes of the (i-dimensional simplices to which the sub-simplex r^-2 belongs. 
In the same way we have written the curvature density R as S nd _ 2 /V nd _ 2 (d) viewing 
it formally as representing some average value in the volume V nd _ 2 (d). With such 
an interpretation one can of course write 

2 



J d d ^R 2 ~ £ V nd _ 2 (d) 



n d -2 



^n d _ 2 Vn d _ 2 (d-2) 
Vn d .M 



(2.24) 



This definition must be interpreted with some care if we want convergence to the 



continuum value for a smooth manifold by successive subdivision ||19|| . We do not 
have to worry too much about these subtleties here since our task in the functional 
integral is not to approximate a given smooth manifold but to select some class of 
manifolds which can be used in an (approximate) evaluation of the integral. From 
this point of view we will use the R 2 term as representing typical higher derivative 
terms which one would have to insert in order to stabilize the Euclidean path integral 
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as explained in the introduction. As is well known from for instance lattice gauge 
theories discretized versions of higher derivative terms are by no means universal. 
It is clear that this point of view is not as beautiful as the original geometrical way 
that Regge viewed piecewise flat manifolds, but our perspective is that a term like 
( |2.24[ ) will probe a universality class of theories which have an effective expansion 
in terms of higher derivative actions like (|2.24| ). 

In the case where we consider the piecewise flat manifolds which can be obtained 
by the process of dynamical triangulation as described above ( |2.24j ) simplifies and 
we get 



d d £y/gR 2 ~ £ o(n d _ 



n d -2 



c d ~ o(n d -2) 
o(n d _ 2 ) 



[2.25) 



This formula has a slight problem with the continuum interpretation since flat space 
does not have a regular tessellation except for d = 2. This means that the term can 
never scale to zero: 

2 



J2 °i n d-2) 



rtd-2 



Cd ~ o(n d -2) 
o{n d _ 2 ) 



> const. N d . 



(2.26) 



If we introduce a scaling parameter a, which is to be identified with the link length, 
and which is going to be scaled to zero, the physical volume V = / d A ^^/g being kept 
fixed, and if we assume canonical scaling of the terms involved of both sides of ( |2.26| ) 
we get (restricting ourselves to four dimensions which have our main interest): 



d^VaR 2 



> 



const. 



DT 



DT 



[2.27) 



This means that the leading term on the lhs Q2.26 ) is just a cosmological constant 
term and by expanding the lhs we see that it also contains an Einstein-Hilbert term 
etc.. Under the assumption of naive scaling we have to write instead of ( |2.25| ) 

2 



Yl °( n d-2) 
nd-2 



Cd - o\n d -2) 



o(n d - 



Co 



+ %R + c 2 R 2 



(2.28) 



Our lattice "i? 2 "-term is thus to be considered as a generalized higher derivative term 
which, when added to the lattice Einstein-Hilbert term, in addition will lead to a 
redefinition of the bare cosmological coupling constant and the bare gravitational 
coupling constant. 

One interesting aspect of the dynamical triangulation approach is that for a finite 
lattice volume it automatically provides a cut off for the Einstein action. This is 
not the case for the conventional Regge calculus where the action can go to infinity 
without the volume diverging. The reason is that the volume (for d > 2) of the 
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d — 2-dimensional sub-simplices can diverge without the corresponding volume of 
the (i-dimensional simplices going to infinity. In the case of dynamical triangulations 
we have (restricting again ourselves to four dimensions) 

- const. ■ N 4 < XX c 4 ~ °( n 2)) < const. ■ N 4 . (2.29) 



//we assume a conventional scaling in the tentative continuum limit we can rewrite 

o 

r . rnnQf r 

(2.30) 



as 



dt a 2 ' 



DT 

where again we have have introduced the link length a, which is going to be scaled 
to zero while the physical volume V = J d 4 £ Jg is kept fixed. 

2.4 Observables 



Let us here discuss the observables (see also |20 1 for a more general discussion) . Due 
to diffeomorphism invariance and the fact that we in quantum gravity have to inte- 
grate over all Riemannian manifolds the observables which are most readily available 
are averages of invariant local operators like the curvature R(x) and suitable con- 
tractions of powers of the curvature tensor, like R 2 or R^uXcn an d the fluctuation of 
these averages. In addition one can discuss and measure so-called fractal properties 
of space-time and finally with some effort define the concept of correlators of local 
invariant operators. 

The simplest observable is the average curvature. If we consider the discretized 
partition function we have 

/ d^^(0R(0 oc -J— £ o(n 2 f-^^A = 2±N 2 - N4 ( 2 .31) 
J v d(d + 1) ^ o(n 2 ) 10 

and since the volume is 

J '"'^ K 4dTT)5 o(n!) = Nt (2 - 32) 

we can define the average curvature per volume as 



In (p.33|) (R) is defined for a single manifold. We get of course the quantum version 
by calculating the functional average of (R) over all Riemannian manifolds, weighted 
by e~ s . This is what we will do numerically. The average curvature is a bulk quantity 
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which will allow us to get a quick survey of the phase diagram of four-dimensional 
quantum gravity. In a similar way we can define (with the drawbacks described in 
the last subsection) (R 2 ) by 



2 = Jj^VjgWQ E n2 o(n 2 ) [(c 4 - o(n 2 ))/o(n 2 )] 2 



Again this average is defined over a single manifold and we have eventually to take 
the weighted average over all manifolds in our ensemble. A quantity which will have 
our interest will be (R 2 ) — (R) 2 . 

A more refined, but related observable is the integrated curvature-curvature cor- 
relation. In a continuum formulation it would be 

X («a) = ( / d% d% y/tffrjRfa) y/g&jRfo) ) " (/ d d i^M)R{i)) 2 . (2.35) 

In a lattice regularized theory one would expect that away from the critical points 
short range fluctuations will prevail, while approaching the critical point long range 
fluctuation might be important and would result in an increase in x{ K 2)- The ob- 
servable x{ K 2) is the second derivative of the free energy F — — In Z with respect to 
the gravitational coupling constant In the case where the volume N4 is kept 

fixed we see that 



X (k 2 , iV 4 ) ~ (N 2 ) N4 - (N 2 ) 2 N4 = - d lnZ }^ N4) (2.36) 



From the above discussion we have to look for points along the critical line &D where 
x(/t 2 , N^j/N^ diverges in the infinite volume limit iV 4 — > 00. 

Another observable is the Hausdorff dimension. One can define the Hausdorff 
dimension in a number of ways, which are not necessarily equivalent. Here we will 
simply measure the average volume V(r) contained within a radius r from a given 
point. In pOf the concept of a cosmological Hausdorff dimension d c h was defined. 



It essentially denotes the power which relates the average radius of the ensemble of 
universes of a fixed volume to this volume: 

(Radius)^ ~ Nl /dch . (2.37) 

From the distribution V(r) we can try to extract d c h- If for large r we have the 
behaviour 

V(r) ~ r dh (2.38) 

we can identify d c h and dh- 
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By the use of Regge calculus it is straightforward to convert these continuum 
formulas to our piecewise flat manifolds. Between two points in the piecewise flat 
manifold there is a geodesic which is a piecewise linear path. Rather than using 
this definition we will use approximations which are much more convenient from a 
numerical point of view and which one expects should be sufficient for the purpose of 
extracting general scaling behaviour and fractal properties: We define the distance 
between two vertices as the shortest path along links which connects the two vertices, 
i.e. it is essentially the number of links of the shortest path since all links have the 
same length. We call this distance the "ni" -distance between vertices and denote 
its value by d\. In the same way we can define a shortest path between two 4- 
simplices as the shortest path obtained by moving between the centers of neighbour 
4-simplices. We call this distance the "n 4 " -distance between 4-simplices and denote 
its value by dj. The dual graph to a given triangulation, obtained by connecting 
the centers of neighbour 4-simplices, will be a 5 -graph and the ^-distance on the 
triangulation will be the ni-distance on the dual 5 -graph. A priori these distances 
are not related and it is easy to find triangulations where they differ vastly for 
specific choices of vertices and associated 4-simplices. But for averages over vertices 
and over triangulations one would expect that they carry the same information 
about the geometry and we shall see that this is indeed the case. 

With the ni- and ^-definitions of geodesic distances on the triangulations it 
becomes trivial to measure numerically relations like (|2.38|) . As we shall see it is 
less trivial to extract in a reliable way the Hausdorff dimension dh- 

Let us finally discuss the measurements of correlation functions in quantum 
gravity. If we have local invariant operators 0\(x) and O^x) (like for instance 
R{x)) we can define a correlation function of geodesic distance d by 

G(d) = J J d%^^)d%^)0 1 ^ 1 )0 2 (^)5(d^ 1 ^ 2 ) - d). (2.39) 

In this definition <i(£i,£ 2 ) denotes the geodesic distance between £1 and £2 for a 
given manifold, but (|2.39|) makes sense also as a functional average, and therefore 
in principle in quantum gravity. When we discuss numerical data we will always 
have in mind the functional average. It might be convenient to divide by a volume 
element to get the dimension of a point-point correlation function. If V(r) denotes 
the volume inside a "ball" of radius r we can write 

dV(r) = V'(r)dr 

and we define 

G(d) = G{d)/V\d). (2.40) 
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In the case where we have a finite Hausdorff dimension the exponential fall off of 
G(d) and G(d) are identical, but it needs not be the case if we have an infinite 
Hausdorff dimension. 



3 The numerical method 

Unfortunately the analytic methods of two dimensional simplicial gravity have not 
yet been extended to higher dimensions. The numerical method of "grand canonical" 



Monte Carlo simulation, which is well tested in two dimensions pl| , |22| , [23], has 



recently been applied to three dimensions [24, 25, 26, 27, 28] and in four dimensions 



|31 , |32] , |33j |34j[ A necessary ingredient for Monte Carlo simulations in simplicial 
quantum gravity is a set of so-called "moves" , i.e. local changes of the triangulations, 
which are ergodic in the class of triangulations we consider. A general set of moves 
in any dimension has been known for a long time ||35| . They are however not optimal 
for numerical simulations. A more convenient set of moves for higher dimensional 



gravity was suggested in |25| . The ergodicity of these moves in three dimensions was 



proved in | 2q| , and the generalization of this proof to d = 4 was given in [p9 j. In d 



dimensions there are d + 1 moves. Their general description is as follows: remove 
an z-dimensional simplex of order d + 1 — i and the higher dimensional simplices of 
which it is part, and replace it by a d — z-dimensional simplex ("orthogonal" to the 
removed i-dimensional simplex) plus the appropriate higher dimensional simplices 
such that we still have a triangulation. Such moves will be allowed provided they 
do not violate ( |2.5| ) and provided they do not create simplices already present, i.e 
which already have the same vertices as the newly created simplices. 

Let us consider four dimensions where there are five moves. The first move con- 
sists of removing a four- dimensional simplex n^old) and inserting a vertex no (new) 
in the void interior and adding links (and the induced higher dimensional simplices) 
which connect no (new) to the five vertices of n^old). In this way n^old) is replaced 
by five new n4's and the total change of AVs is 

AAo = 1, AiVi = 5, AN 2 = 10, AA" 3 = 10, AA" 4 = 4. (3.1) 

The second move consists of removing a three-dimensional simplex n 3 and the two 
n4's sharing it, and then inserting a link ("orthogonal" to n 3 ) and associated n 2 's, 
ns's and n^s. The total change of iVj's is in this case 

AAo = 0, AiVi = 1, AA" 2 = 4, AA" 3 = 5, AiV 4 = 2. (3.2) 

The third move is a "self-dual" move where AiVj = 0. It consist of removing a 
triangle of order three and associated higher dimensional simplices and inserting the 
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"orthogonal" triangle and its associated higher dimensional simplices such that we 
still have a triangulation. The fourth move is the inverse of the second move, while 
the fifth move is the inverse of the first. 

The change in the action induced by these moves can now readily be calculated 
and we are in a position to use the standard Metropolis updating procedure. The 
weights required for detailed balance are easily determined. Let us only remark 
here that the nature of the problem naturally suggests to use indirect addressing 
by pointers since there is no rigid lattice structure. In addition we found it most 
efficient to keep pointers to vertices of order five, links of order four and triangles 
of order three since these are the ones used in the updating. Since programs of the 
above nature are not well suited for vectorization it is optimal to run them on fast 
workstations. 

Since we are forced to use a grand canonical updating where the volume of the 
universe N 4 is changing, it is convenient to use the technique first introduced in 
30] and used successfully in the simulations in three-dimensional gravity |f24, 27] . 



It allows us to get as close to a canonical simulation as possible and it provides 
at the same time an estimate of the critical point k c 4 {k2) for a given value of the 
coupling constant k^. The idea is the following: Assume we want to perform a 
measurement at some fixed value N 4 (F). The task is to constrain the fluctuations 
of N 4 to the neighbourhood of N 4 (F) without violating ergodicity. First we make 
an approximate estimate of the critical point k 4 , which we denote k 4 (N 4 (F)). It 
can in principle depend on N 4 (F). Next we choose the the actual n 4 used in the 
simulation as a function of the value of N 4 in the following way: 

K ( N ) -I <{N 4 {F))-Ak 4 for N 4 <N 4 (F) 

k 4 {1\ 4 ) - I K c( N4 ( F ^ + A/?4 for Ni > Ni ( F y V-V 

For sufficiently large values of N 4 (F) and small values of Ak 4 we will get an expo- 
nential distribution of N 4 s peaked at N 4 (F): 

By monitoring Ak 4 we can effectively control the width of the distribution of N 4 
around N 4 (F) without violating the principle of ergodicity. A measurement of the 
exponential distribution also allows us to determine k\. If the exponentially fall 
off is different above and below N 4 (F) it means that k c 4 is different from our guess 
k 4 {N 4 (F)) and we can use the optimal value in the next run. Further the mea- 
surements of k% for different values of N 4 (F) allow us to extract the subleading 



17 



correction to the distribution. Assume that the partition function has the form: 

Z(k 2 , k 4 ) = E ^(«2» ^ 4 ) e - K4jV4 (3.5) 

where 

Z(k 2 , iV 4 ) ~ jv/^V^ 2 ^ 4 (1 + 0(1/N 4 )) . (3.6) 

By our method the critical point k c 4 determined by measurements in the neighbour- 
hood of various iV 4 (F)'s would lead to: 

^(measurement) = k° 4 {k2) + f^. . (3.7) 

and a determination of the entropy exponent 7(^2), which in fact governs the volume 
fluctuations of the system. 



4 Numerical results 

The measurements of average curvature etc. were performed for different values 
of iV 4 : 4000, 9000, 16000 and 32000. The number of attempted updatings were 
of the order 10, 000 x N 4 (sometimes considerably longer at critical points where 
thermalization was slow). The phase diagram was scanned varying the (inverse) 
bare gravitational coupling constant n 2 and the U R 2 " coupling constant which we 
denote h. For each value of K2 and h and each value of N 4 this implies a fine- 
tuning of the value of the bare "cosmological" coupling constant k 4 to its critical 
value k 4 (k2, h, JV4). As explained in the last section it is convenient to perform the 
simulations in the neighbourhood of some fixed value of the volume so we choose a 
specific N 4 (F) and limit the fluctuations in volume to some neighbourhood of N 4 (F). 
If we decide to perform the measurements after a given number, n, of Monte Carlo 
sweeps, in practise we perform the actual measurement the first time, after the n'th 
sweep, the system passes a state where the value of N 4 is equal to N 4 (F). This is 
what we mean when we say that the measurements were performed for a given value 
oiN 4 . 

The result for the average curvature is shown in fig. 1 for h = over a large range 
of K2, while fig. 2 shows the average curvature for different values of h, ranging from 
h = to h = 20 and N 4 = 16000. It is not possible with our choice of R 2 term to 
increase h further since the acceptance rates in the Metropolis updating become too 
small. 

We observe the following: In case there is no coupling constant except the cos- 
mological coupling constant the average curvature (R) is zero, which is a nice result 
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since it shows that the selection of manifolds in the context of dynamical triangula- 
tions has no bias towards positive or negative curvature. In case we take k 2 positive 
the average curvature will be positive (this corresponds to the conventional sign of 
the gravitational coupling constant). If we take k 2 negative ("anti-gravity") the 
average curvature will be negative. 

For h = and k 2 ~ 1.1 we see a change towards large positive curvature. The 
same change is seen for h > 0, only we have to go to larger values of k 2 when 
h > 10. For a fixed positive value of k 2 the curvature decreases with increasing h as 
expected. However, as discussed above, the limit h — > oo does not really correspond 
to zero local curvature due to the special form of the our discretized "i? 2 "-term. The 
absolute minimum value of the discretized term corresponds to a constant negative 
curvature: R = —0.046. 

In fig. 3 we have shown the susceptibility defined by (|2.36| ). It can be measured 
directly, or as the derivative of the average curvature. We have used the second 
method, but have also measured the susceptibility directly, with comparable results. 
One sees a clear peak which grows somewhat with volume. This could be taken as 
a sign that the system becomes critical in this region, although the system size is 
not big enough to exclude the possibility of a phase transition of higher order. 

An independent signal of criticality is found by looking at the observable (R 2 ) — 
(R) 2 . In fig. 4a we show its behaviour as a function of k 2 for h between and 20. For 
small h values we see a clear peak with a volume dependence in the region where the 
k 2 susceptibility has a peak too. We note the clear asymmetry between the two sides 
of the peak, especially for h = (fig. 4b). This explains why the position of the peak 
seems to be shifted towards smaller values of n 2 when compared to the susceptibility 
curve. Here again our system is not big enough to exclude the possibility that the 
increase with volume is only a finite size effect and that eventually we shall observe 
only a discontinuity in the derivative at a critical point, which again could signal a 
phase transition of higher (perhaps 3rd) order. The signal deteriorates somewhat for 
large values of h, contrary to the susceptibility signal. We can draw a critical line in 
the (k 2 , /i) -coupling constant plane, fig. 5, and fig. 6 shows the average curvature at 
the transition point as a function of h for = 4000, 9000 and 16000. We see that 
the average curvature gets smaller when h increases and also when N4 increases, but 
for all the values we have been able to probe we have 



and we can not conclude that (R) c > is a finite volume effect. The data indicate 
rather that (R) c remains positive even for h — > 00. For h > 12 there seems to be a 



(R) c = (R(K c 2 ,h)) >0 
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qualitative change in the distributions, which might indicate a different transition, 
but we have not found that it cured the problems of the h = situation. 

The computer simulations in three dimensions revealed a similar situation: A 
transition and a (R) c > 0. In three dimensions there was a very strong hysteresis 
in the same transition, favouring a first order transition. Here we have not seen 
the same strong hysteresis. For N± = 4000 there was no problem moving from one 
phase to the other. For larger N4 we have observed very slow thermalization and 
huge fluctuations in geometry close to k>2 = ^2c, but it did not present itself as clear 
hysteresis. 

Let us now explore the change in geometry along the critical line. Above we 
have defined the "geodesic" link distance d± between vertices and the four-simplex 
distance between four-simplices. The average values describe typical radii of 
our universes. They are shown in fig. 7 {(di)) and fig. 8 {(d^)). Although the d± 
distances are approximately six times larger that the di distances they clearly behave 
qualitatively in the same way and reveal a drastic change in the geometry as we pass 
the critical region of K2. The nature of the change seems to be independent of h. 

The typical universes generated by the computer simulations have small radii, 
almost independent of the volume if we are below the critical k 2 region. After we 
have passed the critical region the radii become quite large and show a very clear 
volume dependence. In fact it seems as if the radius grows almost linearly with 
volume. Qualitatively this implies that the Hausdorff dimension is large below the 
critical region and small (in fact close to one) above the critical region. We have 
not attempted to determine the larger Hausdorff dimension. The growth in radius 
with volume is so small that one has to go to much larger volumes in order to do it 
in a reliable way.f] 

We get a nice representation of the change in geometry between the two phases 
by showing the actual distribution of geodesic length in the universes. This is done 
in fig. 9 for the distances for h = 0. If space-time has a fractal structure with 
some Hausdorff dimension dh the distribution should be like 



In fig. 9 we have shown four curves which correspond to K2 = 0.9, 1.0, 1.1 and 1.2 

in the critical region. Since the transition is smooth and extrapolates from large to 

3 It has been argued that one should not use the "geometrical method" advocated here as a 
measure of the Hausdorff dimension, but extract it from the random walk representation of the 
massive propagator, since this seems to give more "reasonable" values. We disagree with this 
point of view. The two methods are mathematically equivalent p6[ and disagreement reflects in 
our opinion the fact that it is not possible to extract the Hausdorff dimensions with the desired 



P(d) ~ d^ 1 




precision. 
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small dh it is of course possible to find a k 2 in the transition region where we get a 
curve quite similar to ( |4.2p with 4 ~ 4 which is of course amusing, but we do not 
consider the value as especially well determined by the numerical simulations. 

Finally the difference in geometry between the highly connected phase for k 2 
below the critical region and very extended phase above the critical region is reflected 
in the curvature distribution. In fig. 10 we have shown the distribution of the average 
curvature per simplex, defined as an average over all ten triangles forming a 4d 
simplex. In addition we plot the values obtained by blocking the value of curvature 
over larger and larger regions (3 and 5) in the d^ distance. The value of h is 
and the k 2 values are chosen in the neighborhood of the phase transition. The 
same qualitative behaviour is observed for larger values of h. In all the cases the 
distribution seems to consist of two parts: one sharply peaked at small curvature 
values and the second rather broad, shifted towards the positive curvature. In the 
highly connected phase the broad part disappears after blocking and already after 
one or two steps the distribution approaches a 5-function. On larger scales we simply 
have a space with constant curvature R. In the other phase it is the peaked part, 
which disappears after blocking and it seems as if the distribution approaches some 
non-trivial limit. 

One disturbing aspect of our results if we want to give the above mentioned phase 
transition a continuum interpretation is the fact that (-R) c > 0. If we take this result 
literally it is difficult to attribute any sensible naive continuum scaling to the system. 
If we introduce a scaling parameter a, which conveniently can be identified with the 
link length in the triangulation the simplest scaling would be one where a — > 
while the volume a 4 A^ 4 was kept fixed. One would then expect the following relation 
between the "bare" lattice curvature and the continuum curvature: 

(R) (lattice) = R(continuum)a 2 (4-3) 

which shows that if R(continuum) should remain finite in the scaling limit a — > 
(R) must scale to zero. It does not. We have found no way to repair this and 
although it is possible to find a scaling 

\(R(K 2 ,h))-R c (K c 2 ,h)\ ~ \k 2 -nit' 1 (4.4) 

its significance is not clear to us due to R c > 0. One possible explanation is that 
our formalism does not admit a tessellation of flat space and that this in some way 
reflects itself in an expectation value of R ? Another possibility is, as mentioned 
above, that the observed phase transition is in fact higher order and the naive scaling 
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relations need not hold although it is very difficult to imagine an alternative scenario 
which would correspond to a physical scaling. 

The appearance of R c > for h = was one of the motivations to look at theories 
with higher curvature terms. The other motivation was to investigate the question 
of universality in the spirit outlined by Weinberg, as explained in the introduction. 
However, our results are negative in the sense that even if R c indeed decreases with 
increasing h, it does not go to zero and for not too large values of h we clearly are in 
the same universality class as for h = 0. Distributions, the nature of the transition 
etc. seem to be the same except for a displacement in and K2. As mentioned 
before, for h > 12 we see a qualitative change in the distributions, which might 
indicate a different transition. This point requires however further investigation. 



5 Discussion 



The notation of a "hot" and a "cold" phase in quantum gravity in d — 3 was 



introduced in p7| . In the hot phase the large entropy of "quantum" universes 
was dominant. These quantum universes were characterized by a large Hausdorff 
dimension and a high connectivity and the hot phase was continuously connected to 
"anti-gravity", where the (bare) gravitational coupling constant is negative. In the 
cold phase extended structures dominated. In fact the Hausdorff dimension seemed 
close to one, suggesting some kind of linear structure. This phase was interpreted 
as representing the dominance of the conformal mode. In the regularized theory the 
action is not unbounded from below, but instead some lattice configurations which 
are pure artifacts without any connection to the continuum will dominate. These 
were the extended structures observed in three dimensions. The interesting question 
was asked, whether it was possible at the transition point to have truly extended 
structures, relevant for continuum physics. In 23] it was shown that the transition 
in three dimensions was of first order, and a continuum limit was ruled out from 
this point of view. 

From a superficial point of view the situation does not look so different in four 
dimensions. We have two phases which we again can call hot and cold. The hot phase 
is continuously connected to the "anti-gravity" region where the bare gravitational 
coupling constant is negative. As in three dimensions the cold phase is characterized 
by an almost linear, extended structure, while the hot phase has a larger Hausdorff 
dimension and much larger connectivity. In the hot phase the average order of 
vertices is much larger and the average curvature changes from being large positive 
in the cold phase to small positive or even negative in the hot phase. However, the 
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nature of the transition seems different in four dimensions. We have not seen any 
true hysteresis, but there are very long thermalization times in the cold phase where 
the linear structures developed. This is in contrast to three dimensions where a 



very pronounced hysteresis was observed p8| . Our data are not incompatible with 
a second order transition, and this opens for the possibility that a continuum limit 
can be associated with the transition. The scenario is from this point of view quite 
nice: In three dimensions a first order transition rules out a continuum limit, but we 
do not really want the continuum limit in a usual sense in three dimensions since we 
would be confronted with the embarrassing question of a three dimensional graviton. 
The physical Hilbert space of pure three-dimensional quantum gravity is most likely 
finite dimensional Q and does not allow for true dynamical fields. The situation in 
four dimensions is probably very different and it is interesting that the discretized 
model seems to hint at such a difference. 

It is also encouraging that our data share some similarities with the results 
obtained by Regge calculus (|| and references therein). As explained in the intro- 
duction the philosophy of the two methods is quite different and it would be a strong 
argument in favour of universality if one manages to obtain the same results by the 
two methods. In Regge calculus one also observes the two phases, and the phase 
with large positive curvature is characterized by very singular spiky configurations. 
They seem similar to our linear structures, which however cannot arise by single 
points moving away from the rest, as is the case in the Regge formalism. The "hot" 
phase is in the Regge formalism characterized by a small negative curvature which 
however (contrary to our results) scales to zero at the critical point. The latest 
results indicate that one actually has a first order transition for h = and only 
for a finite h > it changes to a second order transition. This transition for finite 
h might have some similarity with the change we have seen for large h, but we still 
have R c > 0. We postpone the discussion of this point to later publication. 

At this point we should emphasize again that we see one major obstacle to taking 
a continuum limit at the critical point and that is the fact that the average curvature 
does not scale to zero. As is seen from the fig.l we have (R) ~ for N4 large in 
the case where there is no gravitational coupling constant. This is a nice result 
as it implies that the measure J2t~s 4 selects positive and negative curvature with 
the same weight, as already noticed. As soon as we add a gravitational coupling 
constant we get an expectation value (R) = Rq(k2), which is essentially a linear 
function of k 2 for k 2 not too large. Since the only interesting critical behaviour takes 
place for k 2 > we need at least a reinterpretation of the scaling limit in order to 



The situation might be different if wc include matter fields. 
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be able to claim we can make contact with continuum physics. One possibility is 
that the expectation value of R c has its root in the missing tessellation of flat space, 
but we feel that would be a surprise in a quantum theory of gravity. A more radical 
point of view would be that there is no four-dimensional quantum theory of gravity. 
Maybe quantum gravity needs matter, as the quantum theory of matter might need 
gravity ? 
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Figure Captions 

Fig. 1 The average curvature (R)(k 2 ) for h = and iV 4 = 4000(0), ^4 = 9000(A), 
iV 4 = 16000(+) and iV 4 = 32000(x). 

Fig. 2 The average curvature (i2)(« 2 ) for h = 0, 10 and 20 and N 4 = 16000. 

Fig. 3 Susceptibility xM for h = and 20 and N A = 4000(0), -^4 = 9000(A) and 
iV 4 = 16000(x). 

Fig. 4a (R 2 )-(R) 2 as a function of k 2 for h = 0, 10 and 20 and N 4 = 4000(0), 9000(A), 
16000(+) and 32000(x). . 

Fig. 4b (R 2 )-(R) 2 as a function of k 2 for h = and iV 4 = 4000(0), 9000(A), 16000 (+) 
and 32000(x). 

Fig. 5 Critical line ^(h). 

Fig. 6 Average curvature at the phase transition for iV 4 = 4000(0), 9000(A) an d 
16000(x). 

Fig. 7 Average d x distance for h = and 20 and for N 4 = 4000(0), 9000(A) and 
16000(x). 

Fig. 8 Average d A distance for h = and 20 and for N 4 = 4000(0), 9000(A) and 
16000(x). 

Fig. 9 Distribution of the d 4 distances for h = and k 2 = 0.9, 1.0, 1.1 and 1.2. 

Fig. 10 Distribution of the curvature per simplex and blocked at a d 4 distance 3 and 
5. 
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